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Abstract 

The canonical-basis HFB method provides an efficient way to describe pair¬ 
ing correlations involving the continuum part of the single-particle spectrum in 
coordinate-space representations. It can be applied to super-conducting deformed 
drip-line nuclei as easily as to stable or spherical nuclei. This method is applied 
to a simulation of the approach to the neutron drip line. It turns out that the 
HFB solution has a stronger pairing and a smaller deformation as the Fermi level 
is raised. However, such changes are smooth and finite. No divergences or disconti¬ 
nuities of the radius or other quantities are found in the limit of zero Fermi energy. 
The nuclear density continues to be localized even a little beyond the drip line. 

PACS: 21.10.Pc, 21.30.Fe, 21.60.Jz 


In nuclei near the neutron drip line, the pairing correlation among the neutrons involves 
signihcantly the continuum (positive-energy) part of the Hartree-Fock (HF) single-particle 
states. In principle, there is no difficulty to treat such nuclei with the Hartree-Fock- 
Bogoliubov (HFB) method, which is the framework to incorporate the pairing correlation 
into mean-held approximations. Indeed, there is no practical problem concerning spher¬ 
ical nuclei niEi. However, deformed nuclei are not so easy to treat. The difficulty 
originates in the huge number of quasiparticle states, most of which are spatially dislocal- 
ized continuum-spectrum states. In the quasi-particle formalism, two methods have been 
used to overcome the difficulty, one using transformed oscillator basis jH] and the other 
using a coordinate mesh but only for axially symmetric nuclei |3]. 

Mathematically, HFB ground states can be expressed in the form of the BCS varia¬ 
tional function. The single-particle states in this expression are localized. They are called 
the HFB canonical basis. This localization makes the level density of the canonical basis 
by far smaller than that of the quasiparticle states because the former is proportional 
to the volume of the nucleus while the latter to the volume of the cavity to discretize 
the positive energy orbitals. The canonical-basis HFB method enables one to obtain 
the canonical orbitals without knowing anything about the huge number of quasiparticle 
states. It can be applied to deformed neutron-rich nuclei without difficulties. 
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The canonical-basis HFB method was originally introdnced for spherical nnclei in 
Ref. jH]. I improved the method and implemented it for deformed nnclei PEI. I also 
fonnd the necessity of momentnm dependence for the contact pairing interactions if one 
employs completely coordinate-space representations (like three-dimensional Cartesian 
mesh, nnlike the radial mesh). 

In qnasiparticle HFB method, the canonical orbitals are obtained from the one-body 
density matrix and thns people have not noticed the existence of a more direct relation to 
the Hamiltonian. The canonical-basis formalism discloses this relation. Namely, canonical 
orbitals above the Fermi level are ronghly the bonnd eigenstates of the pairing Hamil¬ 
tonian. It is not the HF Hamiltonian which generates them. This hnding is helpfnl to 
nnderstand the shell structnre in the continnnm part of the spectrnm [ 7 ]. 

Now, let me show a resnlt of a calcnlation performed with the canonical-basis method. 
It is a simnlation of the approach to the nentron drip line. The system is the N = Z = lA 
nnclens. Instead of increasing the difference N — Z,l modify the parameters of the mean- 
field interaction. Namely, to of the Skyrme force is increased (toward zero from below) 
to raise the Fermi level while t 3 (> 0) is decreased so that the satnration density of the 
symmetric nnclear matter is nnchanged. This approach is taken only becanse the present 
version of my compnter program is designed for N = Z systems. I will examine the 
adeqnacy of this approach in fntnre. 

The interaction in the mean-held channel is the Skyrme SHI force jH] withont the 
spin-orbit term. The conlomb force is also tnrned off. Owing to the omission of these 
two interactions, the single-particle states are fonr-fold degenerated. I take into acconnt 
70 canonical basis states in each of the fonr spin-isospin sectors. The parameters of the 
pairing interaction El are = —880 MeV fm^, Aic = 2 fm and pc = 0.32 fm and 
Pc = oo. 

A three-dimensional Cartesian mesh representation is employed to express the single¬ 
particle wavefnnctions withont assnming any spatial symmetries. The mesh spacing is 0.8 
fm while the edge of the cnbic cavity is 40 fm. 

Fig. □ shows how the HFB gronnd state changes as the Fermi level (A) rises. The 
pairing gap is enhanced almost by a factor of two at the drip line (where A = 0 MeV) 
compared with the solntion for the original SHI force (A = —11.7 MeV). The qnadrnpole 
deformation j3 is decreased by the enhanced pairing. The nnclens has a large prolate 
deformation at A = —11.7 MeV bnt becomes spherical for A > —5.5 MeV. On the other 
hand, the r.m.s. radius does not change so much. Its increase is only 35 % even at the 
drip line. Dislocalization of the density occurs not at A = 0 MeV but at higher A (0.8 
MeV). One can see only smooth changes at A ~ 0 MeV. 

Fig. El shows the energies (expectation value of the HF Hamiltonian) of canonical- 
basis states. One can see that discrete bound states are obtained for both negative and 
positive energies. For A > —5.5 MeV, the nucleus becomes spherical and the levels are 
degenerated. At A ~ —3 MeV, the orbitals are (from the bottom) s, p, s, d, p, f, s, d. 
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Figure 1; The average pairing gap (top), the quadrupole deformation parameter f3 (mid¬ 
dle), and the root-mean-square radius (bottom) plotted versus the Fermi level. 
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g, etc. Here again, there seems to be no violent changes at A ~ 0 MeV. Positive-energy 
localized s orbitals begin to spread over the cavity only for A > 300 keV. 

Considering general properties of HFB solutions jTj, the true ground state must be a 
dislocalized state. The reason for the appearance of a localized solution seems to be as 
follows. The dislocalization of an orbital requires the increase of toward 1 because 
otherwise the orbital is conhned in the pairing potential, which is usually very deep com¬ 
pared with the size of the kinetic energy term of the pairing Hamiltonian. However, since 
v‘^ = 1 corresponds to weaker pairing correlation and larger total energy, a dislocalized so¬ 
lution is not necessarily reached in the iteration process of the gradient method in certain 
circumstances like when the Fermi level is positive but low. In contrast, with methods 
based on the diagonalization of the quasi-particle Hamiltonian, a direct jump to = 1 
can take place and it results in the dislocalization of the corresponding orbital as soon as 
it becomes energetically favorable. 

The localized HFB solutions for positive Fermi energies obtained in the canonical-basis 
formalism may be used as rough approximations to the nuclei just beyond the neutron 
drip line. 
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